Joule expansion of a pure many-body state 
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We derive the Joule expansion of an isolated perfect gas from the principles of quantum mechanics. 
Contrary to most studies of irreversible processes which consider composite systems, the gas many- 
body Hilbert space cannot be factorised into Hilbert spaces corresponding to interesting and ignored 
degrees of freedom. Moreover, the expansion of the gas into the entire accessible volume is obtained 
for pure states. Still, the number particle density is characterised by a chemical potential and a 
temperature. We discuss the special case of a boson gas below the Bose condensation temperature. 
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Irreversible processes can be described within the 
framework of quantum mechanics by taking into account 
a large amount of degrees of freedom. A small number 
of these degrees of freedom are involved in the consid- 
ered process and the other ones are ignored. Most of the 
studies are concerned with composite systems consisting 
of a subsystem of interest coupled to one or several large 
reservoirs. For example, the decoherence and relaxation 
into thermal equilibrium of a quantum system is obtained 
by coupling it to a large heat bath In the context of 
transport, the quantum conductor under study is main- 
tained out-of-equilibrium by connecting it to large free 
particle reservoirs Q • The interesting physics occurs gen- 
erally in the subsystem and the reservoirs degrees of free- 
dom are traced out Q . The reservoirs are assumed large 
enough so that their states remain essentially unchanged 
during the studied process. It seems then natural to sup- 
pose that they are initially at thermal equilibrium. 

It is often argued that the equilibrium states of the 
reservoirs result from their coupling to a thermal super- 
reservoir which is not taken into account explicitly in the 
model. Recently, it has been shown that such an ad hoc 
assumption is actually not necessary to understand the 
relaxation of a quantum system into thermal equilibrium 
@) H, 0]- It nas been found that a boson bath ini- 
tially in a pure state of macroscopically well-defined en- 
ergy can induce this relaxation [7[. The thermalisation 
process has thus been obtained as a consequence of the 
pure quantum-mechanical description of a truly isolated 
composite system. In this approach, the temperature of 
the asymptotic equilibrium state of the interesting sub- 
system is not introduced by hand by thermal averaging 
over the bath initial states. It is determined by the den- 
sity of states of the bath and the macroscopic energy of 
its initial pure state. 

In this Letter, we consider an isolated perfect quan- 
tum gas confined in a finite region of space. The gas 
particles do not interact with each other or with envi- 
ronmental degrees of freedom, they are only subject to 
a static confining potential. We show that the textbook 
example of an irreversible process, the Joule expansion, 



is a direct consequence of the quantum mechanics prin- 
ciples. More precisely, we study the time evolution of 
the number particle density ensuing from an initial pure 
many-body state situated in a subregion of the total ac- 
cessible volume. Following Refs. [ad, 0], we consider 
pure states of macroscopically well-defined energy Con- 
trary to the usually studied case of composite systems, 
the interesting and ignored degrees of freedom are here 
part of the same many-body system. Consequently, the 
duration of the investigated process increases with the 
system size and we must study a finite system. However, 
we obtain, in the thermodynamic limit, a clear expansion 
of the gas into the entire accessible region for times much 
shorter than the Poincare recurrence time of the system. 

To simplify, we first restrict ourselves to the case of a 
one-dimensional perfect gas confined in a box of length L. 
The N indistinguishable particles, bosons or fermions, of 
mass m constituting the gas are described by the Hamil- 
tonian 



Hr = - 



±- [ L ' dx^{x)^ = Y j $- 



c k c k 



(1) 



where ip^{x) and c\ create, respectively, a particle at po- 
sition x and in the single-particle eigenstate k. The sum 
runs over the wavevectors k = mr/L where n is a positive 
integer. We use units in which h = ks = 1- In the fol- 
lowing, we study the particle number density of the gas 
which can be written as 



p(x,t) = Y,e ik *(p k (t)) 



(2) 



where the sum runs over both positive and negative 
wavevectors k and (...) denotes the average with respect 
to the initial gas state \tp). The operators p k are given 
by 



1 \" Jk(2k'-k)t/2m~j x 



2L 



k'^k'-k 



(3) 



k'^O.k 



where Ck = sgn(fc)c|fc|. We remark that p\ = p_ k = p k 
and (po(t)) = N/L for any state \tp). 
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FIG. 1: Particle number density as a function of space for 
e/L = 0.34567 and times t/mL 2 = 0, 5 x 1(T 4 and 0.1 for (a) 
a ID fermion gas at T = 10 5 /mf and exp(-/i/T) = 0.1 (b) 
a 3D boson gas at T = or, equivalently, a single particle. In 
this last case, the first two curves cannot be distinguised. 



We assume that initially the N particles are confined 
in an interval [0,£] where £ < L, and that is a state 
of macroscopically well-defined energy E. We are inter- 
ested in the subsequent evolution of the density @ in the 
thermodynamic limit N 3> 1 with finite mean distance 
IJN and energy E/N . The initial gas state reads 



W) = E i>*\a) 



(4) 



\ot)en E 



where \a) refers to the TV-particle eigenstates of the 
Hamiltonian Hg given by ([T]) with the upper limit L 
replaced by I. The wavevectors of the corresponding 
eigenmodes are denoted by q = niz/1 where n is a pos- 
itive integer. An eigenstate \a) corresponds to a set of 



occupation numbers {n q } obeying 



fcrmions, 



q >o n i 



= N. 



is restricted to the values and 1. 



For 
The 



Hilbert space He is spanned by the states \a) satisfying 
E < J2 q >o n qQ 2 < E + SE where SE is much smaller 
than E but much larger than the maximum level spacing 
of H(. In the thermodynamic limit, the dimension D of 
He is practically proportionnal to SE and a density of 
states n(E,N) = D/SE can be defined. Moreover, this 



density satisfies the Boltzmann's relation 

' E £ 



In [n(E, N)SE] ~ Ns 



N' N 



(5) 



where s — S/N is the gas entropy per particle. As is 
well known, for a perfect gas, the relation ([5|) is a direct 
consequence of the exchange symmetry principle Q . 

We now show that, in the thermodynamic limit, almost 
all normalised states (jlj lead to the same particle number 
density. To obtain this result, we use, following Refs. 
[HjHQj the uniform measure on the unit sphere in TLe 
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More precisely, we will evaluate the Hilbert space av- 
erage and variance of {pk{t)} following from this nor- 
malised distribution. To do so, we note that, in 
the limit N 3> 1, the reduced distribution of a fi- 
nite number q of components ip a n), ■ ■ ■ ,^0(5), equals 
Ilp=i -Dexp(— f|V'a(p)| 2 )/ 7r - With this Gaussian distri- 
bution, we obtain the Hilbert space average 



(&(*)> 



2LD ^ 



E 

: k';q>0 



e- lkk ' t/m {q\k'){q\k'~k)n q 
(7) 

where (q\k) = 2(L£)~ 1 / 2 J dx sin(ga:) sin(fcx) and tk — 
k 2 /2m. This expression simply states that the Hilbert 
space average of the expectation value (pk) is equal to 
the microcanonical average at energy E of the opera- 
tor pk- It can be further simplified using the following 
standard arguments Q. The microcanonical probabil- 
ity distribution of the occupation number n q is P{n q ) = 
n(E - n q e q ,N - n q )/n(E,N) where e q = q 2 /S 
is the density of states of the Hamiltonian Hi 
This density obeys the Boltzmann's relation ([5]) with the 
corresponding entropy s. Expanding this entropy in the 
energy n q e q <C E and number n q <c N and taking into 
account that s = s in the thermodynamic limit, result in 
P(n q ) oc exp[— n q (e q — fi)/T] where T and \i are given by 



e q c q c q- 



L = d E S(E,N), t 



-d N S(E,N). 



(8) 



Consequently, the microcanonical average numbers in ([7]) 
can be replaced by the Bose-Eintein or Fermi-Dirac oc- 
cupation function at the microcanonical temperature T 
and chemical potential p, determined by the intensive pa- 
rameters E/N and t/N of the many-body state ^j. 

We define the Hilbert space variance a 2 of (pk) as the 
average of | {pu) — (pk} \ with respect to the measure ((6]). 
We find 



G 2 = 



- y 

D 2 ^ 



\a).\f))£_n t 



HpMW < ^(Pk(tf)E (9) 



3 




4 8 12 

t (10" 4 ml_ 2 ) 

FIG. 2: Particle current J as a function of time for fermions, 
a fugacity of 0.1, T = 10 5 /ml 2 and I = 0.34567L. The dotted 
and dashed lines are, respectively, the Landauer current and 
the contribution of the wavevector k = n/L. 

where (A)e = J2\ a )en E ( a \A\a)/D denotes the micro- 
canonical average at energy E of the observable A. The 
upper bound is simply obtained by replacing the sum 
over the states \(3) € He by a sum over all the N- 
particle eigenstates \j3) of the Hamiltonian Hi. From 
the above arguments, one finds that the microcanoni- 
cal distribution of two eigenmode occupation numbers 
is P(n q )P(n q ') with T and /i given by 1(5)). As a con- 
sequence, the microcanonical average on the right side 
of the inequality ^ is equal to a finite grand-canonical 
average and hence, as D ~ exp(iV), the variance a 2 van- 
ishes exponentially in the thermodynamic limit. There- 
fore, for almost all states (j4|), the particle number density 
is given by {2]) with the Fourier coefficients 

= IT E (#0(#'-fc)/(<?K ifefe ' t/m (io) 

fc';<j>0 

where f(q) = [exp(q 2 /2mT — p/T) =p l]" 1 depending on 
the bosonic or fermionic nature of the particles. To de- 
scribe a gas state of macroscopic energy E, other choices 
than (j4|) are possible. For example, due to the exponen- 
tial TV-dependence (JSJ) of the density n(E, N), one obtains 
the result (|T0|) for almost all states of the Hilbert space 
spanned by the eigenstates satisfying J2 q n q = N 

and J2 q n q e q < E. We also remark that the above deriva- 
tion is not restricted to the observables pt- The result 
(fl~0|) can be generalised to any n-particle observables for 
n -C N. The expression (|10| holds for higher-dimensional 
systems, provided that the confining potential is indepen- 
dent of a; G [0, L]. More precisely, the number of particles 
p{x)dx between x and x + dx evolves according to (fl~0|) 
with / replaced by a sum over transverse eigenmodes. 
We consider the case of macroscopic transverse dimen- 
sions towards the end of the paper. 

We now discuss the time evolution of the density p. 



For that purpose, it is useful, using the properties of the 
Dirac comb function [§], to rewrite (fT0|) as 

(Pk(t)) = -r-r} dx dx'S(x-x'-\ h 2nL) 

AnL w J ~ e J - e m 

xe ikx r)F(x-r}x' + 2p£) (11) 
??=±i 

where F(x) = J dqe tqx f(q) is the Fourier transform of 
the grand canonical occupation function /. The function 
F does not depend explicitly on the length £, it depends 
only on the intensive parameters E/N and £/N via the 
definitions ([8]). It assumes its maximum value 2wN/£ at 
x = and vanishes for large x. As an example, consider 
the classical Maxwell-Boltzmann limit, A <c l/N where 
A = [ttN / Em) 1 / 2 is the de Broglie thermal wavelength 
{E = NT/2 here), where F(x) oc exp(-7ra; 2 /A 2 ). The 
expression (fTTj) clearly shows that (pk) is periodic with 
period tk — 4Lm/k and (pk(tk/2)) — ±(pk(0)) depend- 
ing on the parity of kL/ir. Consequently, the density 
p is periodic with period AmL 2 /it and p(x, 2raL 2 /n) = 
p(L — x,0). Times of the order of mL 2 are practically 
inaccessible. For example, for He atoms and a length 
L ~ 10 cm, the time period of the particle number den- 
sity is of the order of a few days. 

The gas expansion is described by the relaxation, for 
times t <C mL 2 , of the Fourier coefficients (f>k(t)} corre- 
sponding to wavelengths of the order of the box length 
L. For these times and wavelengths, we obtain 

up to a correction of order L _1 . The macroscopic wave- 
length components of the density p are of order unity 
at initial time and then relax according to (fT2| . The 
gas evolution is hence essentially an expansion into the 
entire accessible volume, as illustrated by Fig[TJ This 
figure clearly shows that this evolution is very different 
from that of a single particle The results displayed 
in Figs [T] and [2] are obtained by numerical evaluation of 
the sum (|10l) . We remark that in the Maxwell-Boltzmann 
limit, the gas expansion described by (fT2")l is characterised 
by an ^-independent time of the order of L(m/fc_BT) 1/ ' 2 
as expected from classical physics dimensional consider- 
ations. In the example already mentionned of He atoms 
in a box of 10 cm length, this time is of the order of 0.1 
ms at T = 300 K and is thus far much shorter than mL 2 . 

Interestingly, the macroscopic wavelength approxima- 
tion (|12p can be interpreted in quasiclassical terms as 
follows. The number particle density resulting from this 
approximation can be written as an integral over p of the 
distribution w(x,p, t) = f(p) ~Y^ n Hi(x—pt/m + 2nL)/2'K 
where Ht(x) = 1 for \x\ < I and otherwise. This 
probability density describes classical particles which 
are initially distributed uniformly in the sub-box [0,^] 
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with momenta p distributed according to the grand- 
canonical occupation function and which evolve freely 
between perfectly elastic collisions with the walls of the 
container [0,L] The gas expansion given by (fT2"|) 

can be characterised by the time evolution of the par- 
ticle current J(t) = —d t J dxp(x,t). For short times, 
J is constant and equal to the Landauer current [l2j 
Jl = de/2n[exp(e/T) =p 1] flowing through a perfect 
conductor from a reservoir at thermal equilibrium, the 
interval [0, £] here, to an empty reservoir, the interval 
[£, L] here, see Fig [2] For longer times, the density p and 
hence the current J are well approximated by the longest 
wavelength terms of ([2]), see Figs. [Hand [21 Finally, the 
current J essentially vanishes. 

Contrary to the quasiclassical approximation Q12p , the 
quantum number particle density determined by (fT0|) 
does not relax into a steady state. As (pk) is time pe- 
riodic with period tk = 4Lm/£;, the density p fluctuates 
constantly, even in the Maxwell-Boltzmann limit. For the 
sake of clarity, we discuss this limit and a length L > 2£. 
From (fTT]) . we deduce that |(pfc(t))| is even and periodic 
with period ife/2. Moreover, we find, for < t < ifc/4, 
\(pk(t))\ ^ /(fc/2)/4Z, for t/t k < £/2L and elsewhere, 
provided that t/tk, \t/tk — £/2L\ ^> A/L where A is the 
de Broglie thermal wavelength. To show that the micro- 
scopic fluctuations of the density p(x, t) do not dissapear 
in the thermodynamic limit, we consider 



M(t) = dx[ p(x,t) - 



N 

T 



= Lj2\{Pk(t))\ 2 (13) 



which is a global measure of the difference between p 
and the uniform density N/L. The distance M decreases 
from the extensive value M(0) ~ (L/£-l)N 2 /L to finite 
values for t 3> mLA. For these times, all the Fourier 
components \(pk)\ are of the order of L _1 and M can be 
written as an integral over k on the intervals [2m (pL — 
£)/t, 2m(pL + £)/t] where p is an integer. Consequently, 
the function p(x,t) — N/L is non- vanishing in finite size 
regions, see Fig[TJ and changes with a characteristic time 
of the order of mLA. 

We finally discuss the case of a gas confined in a tridi- 
mensionnal box of length L and rectangular cross-section 
S ~ L 2 . We assume that the gas is initially in a state of 
macroscopically well-defined energy E and situated in a 
sub- volume V = S£ between x = and x — £ < L. The 
main difference with the unidimensional system is that 
here, for bosons, the average occupation number No of 
the sub-box V lowest eigenmode is macroscopic for E be- 
low the Bose condensation energy Eb oc N(N/V) 2 ^ 3 /m 
8j. We note that, in this case, the microcanonical prob- 
ability distribution of this occupation number is not of 
the grand-canonical form 13(. Whereas one eigenmode 
can be macroscopically populated here, the result p = ~p 
remains valid since the dimension D in ^ grows ex- 
ponentially with N. For fermions and for bosons with 



E > Eb, the gas expands according to ((TTj) with F re- 
placed by F 3 d(x) = S J dqqsm(qx)f(q)/2nx. For bosons 
with E < Eg, the density p — pbec + Pje is the sum 
of two contributions corresponding, respectively, to the 
condensed and non-condensed bosons. The contribution 
Pje is the zero chemical potential limit of the above case. 
The other contribution is very different since pbec(x, t) 
equals Nq times a single-particle position probability 
density, see FigQ] Therefore, though we consider pure 
states, the behavior of the boson gas is radically differ- 
ent below and above the Bose condensation temperature 
Tb^Eb/N |. 

In conclusion, we have obtained an irreversible evolu- 
tion of physically relevant degrees of freedom of a many- 
body system without any coupling to environmental de- 
grees of freedom. More precisely, we have shown that 
an isolated quantum perfect gas confined in a finite re- 
gion of space tends to fill uniformly the total accessible 
volume, with the noticeable exception of the low energy 
3D Bose gas. Moreover, whereas we consider pure gas 
states, the gas expansion is characterised by a chemical 
potential and a temperature and is well described by the 
Landauer formula for short times. It would be interesting 
to study the influence of interactions between the parti- 
cles and with an environment. For a weak coupling to an 
environment, the behavior of the number particle density 
should be well described by our results for short enough 
times. The particle-particle interactions could induce a 
complete relaxation of the single-particle density matrix 
into a grand-canonical equilibrium state. 

We thank R. Chitra for helpful discussions and a care- 
ful reading of the manuscript. 
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